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We use the cluster variation method (CVM) and Monte Carlo simulations to investigate the phase structure 
of the 3d gonihedric Ising actions defined by Savvidy and Wegner. This model corresponds to the usual three- 
dimensional cubic Ising model with nearest, next to the nearest, and plaquette interactions in the region with 
degenerate lamellar ground states. The picture of the phase diagram given by the CVM is in good agreement 
with the results of Monte Carlo simulations, and it is shown that the gonihedric model is in the same universality 
class of the ordinary three-dimensional Ising model. 



The 3d Ising model with interactions extended 
to next-to-the-nearest neighbouring (NNN) sites 
and 4-spin plaquette interaction has a very rich 
phase diagram relevant for the description of 
physical systems of interacting interfaces [QJ|] . In 
this paper we discuss the critical behaviour of this 
model in a particular region of parameters first 
explicitly considered in || in relation with string 
theory 

Ising variables can describe surface configura- 
tions in a simple way to explain: for any spin 
configuration there is on the dual lattice a set of 
closed interfaces separating the domains of spins 
with different sign. If 0a,PctPi respectively are 
the energy cost for an unit area, two connected 
plaquettes forming a right angle, and four pla- 
quettes intersecting at a common link, the corre- 
sponding Ising model is given by 



-H - Ji E(ii> + h E<(i,i» 



(1) 



with J x = I3 C + (Pa + Pi)/2,J 2 = -(20c + 
/3/)/8, J 4 = (20 c ~ Pi)/8 §. In eq. (@) the three 
sums respectively are on the nearest neighbours, 
the NNN spins and the plaquettes of the cubic 
lattice. The gonihedric model is the model ([l]) 
with p A = (J 2 /Ji = -1/4) and C = 1- The 
choice Pa = means that the amount of sur- 
faces in the system is not controlled by any ex- 



ternal parameter, while Pc — 1 gives a prefer- 
ence to flat surfaces. The parametrization J\ = 
2k, J2 = —k/2, J4 = (1 — k)/2 is often used; the 
parameter k = Pi J A + 1/2 can interpolate be- 
tween completely non-interacting (k = 1/2) and 
self-avoiding surfaces (k = 00). 

The ground state degeneracy is the character- 
istic of the gonihedric Ising model: all the pos- 
sible sequences of "+" and "-" spin planes (for 
J4/J1 > —1/4) have the same minimum energy. 
The situation changes for Pa 7^ 0: a positive Pa 
selects the ferromagnetic ground state (no inter- 
face allowed) while at negative Pa the ground 
state is a sequence of planes with spins of alter- 
nate sign (maximum amount of flat interfaces). 

The finite temperature behaviour of the model 
can be easily studied by means of the mean field 
approximation, which predicts that the ground 
state symmetry is preserved at finite tempera- 
ture (that is, the first order line between lamel- 
lar and ferromagnetic phases in Fig. 1 would 
be strictly vertical, independent of temperature) 
and, at least for k > 3/4, the ferromagnetic and 
lamellar phases are separated from the disordered 
one by two second order lines (both turning first 
order for k < 3/4) which meet in a bicritical point 
at J 2 /Ji = -1/4. 

Recently, several studies based on high tem- 
perature expansions || and Monte Carlo simu- 
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lations |fj] addressed the issue of the universality 
class of the gonihedric model (the existence of a 
phase with long-range order has been analytically 
proved in Q for the case J4 = 0), formulating 
the conjecture that the model might be in the 
Onsager (i.e. two-dimensional Ising) universality 
class. Further, longer simulations ||], however, 
found estimates for critical exponents which were 
seriously incompatible with both the conjecture 
and the previous numerical estimates. 

We have studied the gonihedric model at fi- 
nite temperature by means of the cluster vari- 
ation method (CVM) |Q, which is a power- 
ful generalization of the mean field approxima- 
tion p!l[ . Our results show that the ground- 
state degeneracy along the line J2/J1 = —1/4 
is broken at finite temperatures and the ferro- 
magnetic phase is stable along this line. In 
Fig. 1 we show the phase diagram in the space 
(l/(/3Ji), J 2 /Ji(J 4 = 0) (f3 = (fcsT)- 1 ) as cal- 
culated in the cube approximation of the CVM. 
We see that the ferromagnetic-lamellar transition 
line bends at finite temperatures towards val- 
ues of J2/J1 lower than -1/4; this bending has 
been confirmed by low-temperature expansion re- 
sults [p|Jl^1 and is very important since, with the 
reasonable assumption that the ferromagnetic- 
paramagnetic line, which is of second order, does 
not exhibit non-universal behaviour, we see that 
the gonihedric model must be in the same uni- 
versality class of the ordinary three-dimensional 
Ising model. 

The lamellar to disordered transition line turns 
out to be first order for J 4 — (this being in- 
dipendently confirmed by the numerical simula- 
tions in (l3)), and hence the correct phase dia- 
gram exhibit a critical end point structure, in- 
stead of the bicritical one predicted by the mean 
field approximation. Nevertheless, CVM results 
for the critical amplitude of the Ising transi- 
tion close to J2/J1 = —1/4, analyzed with Pade 
approximants p4| sh ow the existence of strong 
cross-over effects |1(J , probably due to the vicin- 
ity of the lamellar spinodal line, which should 
be responsible for the contrasting results for the 
critical exponents found in the simulations ^J||. 
Other previous simulations [jl3[ were also not able 
to give definitive estimates of critical exponents. 



At J2 / Ji = — 1/4 the critical value of 1 / (fl J\ ) is 
1.171. It can also be shown by CVM calculations 
that for other values of J4 the ferromagnetic- 
paramagnetic transition becomes of first order 
close to the lamellar transition with a tricritical 
point at J2 < 0. In the fc-parametrization the 
tricritical point is observed for k < 0.86. 
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Figure 1. The phase diagram of model (|) at 
J4 = 0. Solid and dashed lines represent first and 
second order transitions respectively. 

In order to understand the source of the dis- 
crepancies between the various simulations, we 
tried also to perform a standard Metropolis 
Monte Carlo on this model. In Fig. 2 we give 
our results for the behaviour of the specific heat. 
Our data have been obtained by averaging over 
2000 decorrelated measures at each value of the 
lattice size L and of the inverse temperature (3. 
Our decorrelation times were 1000 full updates of 
the lattice in the case of the largest value of L. 
These times are substantially longer than those 
reported in |]||, but we found them necessary 
(at least with L = 18) in order to get stable av- 
erages. This suggests that the simulations per- 
formed in Q with much larger lattices are defi- 
nitely too short. 

We fitted our data with the function C(j3) = 
C ma x/[fr(/3 — flc) 2 + 1] and the values of C max and 
Pc obtained for each lattice size are listed in Ta- 
ble 1. A fit of the position of the peak with the 
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function [3 C — /3 cr it + aL~i gives the estimate for 
the critical inverse temperature /3 C rit = 0.4370 ± 
0.0001 corresponding to l/(/3J%) = 1.144 to be 
compared with the CVM value 1.171 and with 
1.136 of the previous simulations ](|. Our esti- 
mate for the exponent l/v is \/v = 1.483±0.001. 
Since the Ising value is l/i/ = 1.594 ± 0.004 || 
we have to go to much larger lattices in order to 
get reliable estimates. 



rameter space by means of the cluster variation 
method, that the three dimensional gonihedric 
model is in the same universality class of the or- 
dinary three dimensional Ising model. Further- 
more, the model seems to exhibit very strong 
crossover effects which pose serious problems to 
numerical simulations. 

We are indebted to Marcia Barbosa, Amos 
Maritan, Michele Caselle, Dino Cosmai and Paolo 
Cea for helpful discussions about this work. 
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Cm ax 
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2.8041 ±0.0237 


0.41182 ±0.00021 


10 


4.3707 ±0.0706 


0.41872 ±0.00018 


12 


6.1137 ±0.0334 


0.42300 ±0.00010 


16 


9.4068 ± 0.0586 


0.42794 ±0.00004 


18 


10.6710 ±0.0883 


0.42937 ±0.00002 



Table 1 

Valus of C m ax and f3 c for different choices of the 
lattice size L obtained by fitting data in Fig. 2. 
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Figure 2. Monte Carlo measurements of specific 
heat in the gonihedric model with k = 1. Circles, 
squares and triangles refer respectively to L = 
12, 16, 18. Solid lines represent the best fit of the 
peaks performed with the function given in the 
main text. 



In summary, we have shown, studying the 
phase diagram of the model in an enlarged pa- 
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